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Abstract 

We numerically study the spatially-resolved NMR around a single vortex in a noncentrosymmetric superconductor 
such as CePtaSi. The nuclear spin-lattice relaxation rate Tj -1 under the influence of the vortex core states is 
calculated for an s + p-wave Cooper pairing state. The result is compared with that for an s-wave pairing state. 
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The heavy fermion superconductor CePtaSi has a 
noncentrosymmetric crystal structure [1]. The lack of 
the inversion symmetry leads to the mixture of Cooper 
pairing channels of different parity, resulting in unusual 
properties observed experimentally [1,2,3]. For exam- 
ple, it is observed that while a Hebel-Slichter coherence 
peak appears in the nuclear spin-lattice relaxation rate 
T^ 1 [2,3,4] suggesting an s-wave like Cooper pairing, 
low-temperature experiments on the London penetra- 
tion depth [3,5] and the thermal conductivity [6] in- 
dicate line nodes in the quasiparticle gap (i.e., uncon- 
ventional superconductivity) . In this context, was 
recently discussed in Refs. [7,8,9] theoretically. 

On the other hand, the spatially-resolved nuclear 
magnetic resonance (NMR) method has been revealed 
to be a powerful experimental technique [10,11]. This 
technique as a probe of the electronic structure with 
spatial resolution is expected to reveal pairing symme- 
try of unconventional superconductors, because in spa- 
tially inhomogeneous systems there appear properties 
specific to the unconventional superconductivity. 

In this paper, we study the spatially-resolved T^ 1 
around a single vortex for an s + p-wave pairing state. 
This pairing state is proposed for CePtaSi in Ref. [12]. 
We consider a system described in Refs. [12,13]. From 
now on, the notations are the same as those in Ref. [8]. 

To obtain Tj -1 in the spatially inhomogeneous situ- 



ation, we calculate the quasiclassical Green function g 
which is a 4 x 4 matrix in Nambu and spin spaces, g 
follows the Eilenberger equation for the noncentrosym- 
metric superconductivity [8,14] 
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where a indicates the strength of the Rashba-type spin- 
orbit coupling, gfc is composed of the antisymmetric 
vector gfe = (— k y , k x , 0) [13], S is the electron spin op- 
erator, and A*, is the superconducting order parame- 
ter. For a static magnetic field in the z direction, the 
expression for \jT\T is given as [8,15] 
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where a — (a v ) and a lj = (a^„), [i,j = {1,2}, fx, v = 
{Ti !}]• The brackets (■ • •) denote the average over the 
Fermi surface. r\ (> 0) is a small constant, which can 
roughly represent the impurity effect. 



Preprint submitted to Elsevier Science 



2 February 2008 



We consider two Fermi surfaces split by the spin- 
orbit coupling agk ■ S, and have the s + p-wave or- 
der parameter >]/ + A sin 6 on one Fermi surface and 
^ — A sin 9 on the other one [12] ($: s-wave, A: p-wave). 
We assume that the difference of the density of states 
and vy between the two Fermi surfaces is small, and 
neglect it here. We use the following empirical form of 
a single vortex, $/T c = aF(t) tanh(r /£(£)) exp(i<f> r ), 
f(t) = £oVi/bF(t), F(t) = tanh(1.74v'lA- 1). Here, 
$ = {*, A}, £ = v F /T c , t = T/T c and the fitting pa- 
rameters a, b. The flux line is along the z axis, r is the 
distance from the vortex center, and <j> r is the angle in 
the x-y plane. We have performed a full numerical com- 
putation as in Ref. [16] and obtained a self-consistent 
vortex solution for the same parameters used in Ref. 
[8]. We then fitted it with the above expression, and 
obtained (o,6) = (0.68,2.2) for * and (1.42,2.2) for 
A. We have also done the same fitting for the conven- 
tional s-wave pairing case, so that (a,b) — (1.76, 2.2). 
We found that the fitting is good for T > 0.1T C , while 
for T < 0.1T C it gives rather smaller core radius than 
actual one. 

We calculate 1/TiT at several positions around the 
vortex, and show the results in Fig. 1 (s + p-wave state) 
and Fig. 2 (s-wave state). Away from the vortex (r — 
10£o ), 1/TiT exhibits the coherence peak just below 
T c in both pairing states, while it is rather smaller in 
the s + p-wave state [8]. With decreasing r up to £o, 
the height of the peak gradually decreases. With fur- 
ther decreasing r below £o inside the vortex core, the 
height of the peak increases again and the value of 
\jT\T at low temperatures increases also. The differ- 
ence between \ jT\T at r = 0.6£o and that at r — 10£o 
in Fig. 1 (s + p-wave state) is similar to the difference 
between higher and lower frequency data of experimen- 
tally observed l/T x T in Refs. [2,3,4]. That is, almost 
same hight of the peak just below T c inside and out- 
side the core, and larger values of 1/TiT at low tem- 
peratures inside the core than those outside the core. 
Finally, we note that 1/TiT increases further and the 
peak shifts to lower temperature side with approaching 
the vortex center (see plots for r — 0.2£o in Figs. 1 and 
2). The difference between s+p- wave state (Fig. 1) and 
s-wave state (Fig. 2) is conspicuous there in the low- 
temperature behavior. It would be due to a difference 
of the vortex core states between these pairing states. 

In conclusion, we calculated spatially-resolved 
1/TiT for s + p-wave and s-wave pairing states. The 
information on 1/TiT deeply inside the core might 
be valuable for experimentally detecting the pairing 
symmetry for a noncentrosymmetric superconductor. 
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Fig. 1. 1/TiT vs T for the s + p-wave pairing state, r is the 
distance from the vortex center. 77 — 0.01T C . 
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Fig. 2. 1/TiT vs T for the s-wave pairing state, r/ = 0.01T C 
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